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G^ I We describe diffractive deeply inelastic scattering in terms of diffractive parton 

distributions. We investigate these distributions in a hamiltonian formulation 
that emphasizes the spacetime picture of diffraction scattering. For hadronic 
O-i. systems with small transverse size, diffraction occurs predominantly at short 

r-| ! distances and the diffractive parton distributions can be studied by perturba- 

tive methods. For realistic, large-size systems we discuss the possibility that 
diffractive parton distributions are controlled essentially by semihard physics 
at a scale of nonperturbative origin of the order of a GeV. We find that this 
possibihty accounts for two important qualitative aspects of the diffractive data 
from HERA: the fiat behavior in (3 and the delay in the fall-off with Q^. 



CERN-TH/99-154 
ETH-TH/99-09 
PSU-TH/207 
May 1999 



1. Introduction 

In hadron-hadron scattering, a substantial fraction of the events are diffractive: one 
or both of the initial hadrons emerges with small transverse momentum in the final state, 
having lost only a small fraction of its energy. Assuming that quantum chromodynamics 
is the theory of the strong interactions, one expects that diffractive scattering is due to 
the exchange of gluons. Since gluons are pointlike objects, the gluon exchange picture 
suggests the possibility of hard diffractive scattering, in which exchanged gluons moving in 
opposite directions participate in a hard process such as jet production, with a transferred- 
momentum scale much larger than Aqcd- Similarly, in lepton-hadron scattering it should 
be possible to have a hard process, deeply inelastic scattering, in which the incoming 
hadron is diffractively scattered. These possibilities were suggested in 1985 by Ingelman and 
Schlein fl]] . Experimental data from both hadron-hadron p|, ^ and lepton-hadron [^ ^ § 
colliders have confirmed the existence of hard diffractive scattering. 

Since hard diffractive scattering contains a hard subprocess, one may ask whether the 
QCD factorization theorem that holds for inclusive hard scattering works also for diffractive 
hard scattering. Investigation of this question indicates that factorization does not hold for 
diffractive hard processes when there are two hadrons in the initial state. In inclusive hard 
scattering there are contributions from particular final states that would break factorization, 
but these contributions cancel because of the unitarity of the scattering matrix when one 
sums over all final states 0. If only diffractive final states are allowed, this cancellation is 
spoiled 1^, 0. However, factorization does hold for diffractive deeply inelastic scattering, 
in which there is only one hadron in the initial state p, 0, [ll[]. It is diffractive deeply 
inelastic scattering that is of concern in this paper. 

In diffractive deeply inelastic scattering, one can measure the diffractive structure func- 
tion dF2^^{x,Q'^,Xp,t)/[dxpdt], where Q^ and x are as usual the photon virtuality and 
the Bjorken variable of deeply inelastic scattering, Xp is the fractional loss of longitudinal 
momentum by the diffracted hadron, and t is the invariant momentum transfer from the 
diffracted hadron. This structure function is often called -^2 • The factorization theorem 
allows us to write 

The function Fa is the hard scattering function, calculable in perturbation theory. It is 
the same hard scattering function as in inclusive deeply inelastic scattering. The function 
dfa^h^/[dxp dt] in Eq. (|1.1|) is the diffractive parton distribution, containing the long distance 
physics. It is interpreted as the probability to find a parton of type a in a hadron of type 
A carrying momentum fraction S, and, at the same time, to find that the hadron appears 
in the final state carrying a fraction 1 — Xp of its longitudinal momentum, having been 
scattered with an invariant momentum transfer t. Both the hard scattering function F and 
the diffractive parton distribution functions df^^^/[dxpdt\ depend on a factorization scale 
11. The 11 dependence of the distribution functions is given by the usual renormalization 



group evolution equation, 

9 In fj? dXp dt a ■^ aa , ^^^ ^^ 

where the kernel Kaa' has a perturbative expansion in powers of as{fi)- The difFractive par- 
ton distributions are defined as certain matrix elements of quark and gluon field operators, 
analogously to the definition of the ordinary (inclusive) parton distributions |T^ . We will 
discuss these definitions later in the paper. 

Our purpose in this paper is to investigate the diffractive parton distribution functions, 
expanding on the analysis reported in Ref. [0. We are interested in the leading behavior 
of these functions when Xp <ti I. In the language of Regge theory, this corresponds to 
looking in the region where the pomeron is dominant over other Regge poles. Although 
the evolution equation for the diffractive parton distribution functions is the same as that 
of the inclusive parton distribution functions, their behavior at a fixed scale /iq that serves 
as the starting point for evolution may be very different from the behavior of the inclu- 
sive functions. The different phenomenology that characterizes diffractive versus inclusive 
deeply inelastic scattering depends entirely on this. 

Of course, the diffractive parton distributions in a proton at the scale /xq are not pertur- 
batively calculable. Notice that the problem lies with the large transverse size of the proton. 
Suppose one had a hadron of a size 1/M that is small compared to 1/Aqcd- Then one 
could compute diffractive parton distributions as a perturbation expansion. In this paper, 
we first consider diffraction of small-size hadronic systems. We study this in detail. Then, 
we discuss how the picture of diffraction changes as we let the size increase. This involves 
nonperturbative dynamics. We explore whether one may extract (at least, qualitative) in- 
formation on the diffractive parton distributions for a large-size system by supplementing 
the computation at a much smaller size scale with a hypothesis on the infrared behavior of 
the diffraction process. 

As a simple case of a small-size hadronic system, we consider a diquark system produced 
by a color-singlet current that couples only to heavy quarks of mass M ^ Aqcd- This 
system gets diffracted and acts as a color source with small radius of order 1/M. In this 
case, the perturbation expansion for the leading 1/xp -^ oo terms in the diffractive parton 
distributions begins at order a^. Although this is a rather high order of perturbation 
theory, we find that the result has quite a simple structure and can be expressed in terms 
of integrals that can be evaluated numerically. 

One can view the problem that we address as being that of diffractive deeply inelastic 
scattering at scale Q from a hadron of size 1/M with Q > M and Xp <^ 1. There are two 
main ingredients in our analysis. 



The first ingredient has already been introduced: the factorization formula (|1 . 1|) . Using 
factorization, we are led to analyze the diffractive parton distribution functions, which 
are simpler than -F^'^. This ingredient is especially important because (as we will see) the 



diffr active gluon distribution makes an important contribution to -F^'^, but this contribution 
is not so easy to analyze systematically without the use of the factorization formula. 

The second ingredient is the physical picture that, in a suitable reference frame, the 
partons that are "measured" in the process are created by the measurement operator outside 
of the hadron and, much later, interact with the hadron. This picture, called the "aligned 
jet" model by Bjorken |]13|, applies when Xp <^1. We will see that it emerges most naturally 
when one works in configuration space using light cone coordinates 

x± = ^4^ . (1.3) 

The calculation becomes particularly transparent when one uses a hamiltonian formulation 
in which the theory is quantized on planes of equal x~ . We will see how this works in 
Sec. 3. 

The plan of the paper is as follows. In Sec. 2 we review the operator definitions for 
diffractive parton distributions and describe their structure at large Xjxp. In Sec. 3 we show 
how to compute these distributions in the light cone hamiltonian formulation of the theory. 
In Sec. 4 we present general properties and numerical results for the distributions. In Sec. 5 
we comment on their evolution and the structure of ultraviolet divergences. In Sec. 6 we 
discuss the relation of the previous calculations with the phenomenology of diffractive 
deeply inelastic scattering at HERA. We give conclusions in Sec. 7. In Appendices A-D we 
give calculational details on certain operator matrix elements, we outline the main steps of 
the covariant formulation alternative to the one presented in the text, and we collect some 
integrals. 



2. Basic definitions and approximations 

In this section, we outline the parts of our analysis that bear on the general structure 
of the result. We begin with the operator definitions of the diffractive parton distributions. 
Then we describe how the diffractive parton distributions break into a convolution of a 
part associated with these operators and a part associated with the wave function of the 
incoming state. 



2.1 Operator definitions 

Let us briefly recall the definition of the diffractive parton distributions in terms of 



matrix elements of bilocal field operators |10|. This is the same definition |14, IH] as 



for inclusive parton distributions except that one requires that the final state include the 
diffractively scattered hadron. Let {j)a-,sa) and (pa'^^a') denote the momentum and spin 
of the incident and the diffracted hadron. Let us adopt the standard notation ^ixp for the 



hadron momentum fraction carried by the parton. For gluons one has 






dxp dt 16 TT^ 27r (3x]pP^ ^ sa " x s , 

X (pa, SA\Fa{Qy^\pA'. SA'] X) {pA' , sa'] X|F,(0, y' , Oy^PA, sa) , (2.1) 

where there is an imphcit sum over j = 1,2 and where Fa{0,y~,O)~^''' is the field strength 
operator modified by multiphcation by an exponential of a line integral of the vector po- 
tential: 

Fa{yy^ = E{yUF,{yy^ , (2.2) 

wherePl 

/ TOO \ 

E{y)=Vexp[-igJ dx~At{y+,x-,y)t,j . (2.3) 

The symbol V denotes path ordering of the exponential. The matrices t^ in Eq. ( |2.3| ) are 
the generators of the adjoint representation of SU(3), {tcjbd = —ifcbd- 



The field operators in Eq. ( |2.1| ) are evaluated at points separated by lightlike distances. 
There are ultraviolet divergences from the operator products. It is understood that these 
are renormalized at the scale /i using the MS prescription. 



For purposes of computations, there is a more useful way ||T4l to write F^K Starting 
with 

FM^^ = EiyU Id^Abiyy - d^ A^iy^ + tgA^+A^iyy (tj J , (2.4) 



we have 



FaivV' = d+iE{y)atAbiyy)-igE{yUA,{yyiQbdA,{yy 

-Eiy)abd^A,iyy + igEiy)abA,iyyA,iyy (Qm ■ (2.5) 



The second and fourth terms cancel. Furthermore, inside the integration in Eq. (|2.1|) , 
d~^ = d/dx~ becomes —il3xpp\\ 

Fa{yy^ ^E{y)ab[-il3xpp\Ab{yy-d^Ab{yy] . (2.6) 

This form is useful because it does not have an A x A term and it does not have a d/dx~ . 
(In the complex conjugate matrix element in Eq. (|2.1|) , the —i(3xpp\ becomes +i(3xpp\.) 

For quarks of type j one has 

dffuiPxp,Xp,t,a) 1 1 1 ^^ r fl+^-ir^ 

dxp dt 16 7r2 27rpXpp\ 2 j^ J ^ £^^ 

X ^^(PA, SAl'^mpA', SA'-, X)7+(pA,, SA'-, X|^(0, y", 0)^^, Sa) , (2.7) 



*Here the sign in the exponent is opposite to that in Refs. [|l2[ ^. The sign choice depends on the 
convention for the sign of g. Here we choose the sign of g so that Z)^ ~ d^+ igA°jta- 



where 

^>,{y) = E{y)mj{y) . (2.8) 

Here E{y) is given by Eq. ( |^.3| ) with, now, the matrices t^. being the generators of the 
fundamental representation of SU(3). 



2.2 General structure at large 1/xp 

As discussed in Sec. 1, the scattering process we consider is initiated by a color-singlet 
current. We specialize the definitions of the previous subsection to the case in which the 
incident hadron A is a special photon that couples to a heavy quark-antiquark pair of mass 
M. This pair couples again to the diffractively scattered photon in the final state. 

In the discussion that follows, we use the same method for the diffractive quark distri- 
bution and the diffractive gluon distribution. For the sake of definiteness, we present the 
discussion in the case of the diffractive gluon distribution. At the end of Sec. 3 we will give 
the extension of the results to the case of the quark. 

We consider the limiting 1/xp -^ oo behavior of graphs for the diffractive gluon dis- 
tribution. We shall find contributions proportional to s"^ at fixed (3. At higher orders in 
perturbation theory, the x~^ would be supplemented by logarithms of Xp. In this paper we 
will limit ourselves to considering the graphs at the lowest order of perturbation theory for 
which x~^ behavior emerges. 

What are the lowest order graphs that can give x~^ behavior? Note that there must 
be at least one gluon in the final state in order to balance the color of the operator that 
measures the gluon distribution. We find that with exactly one gluon in the final state 
there are x~^ contributions at order a\. An example of a contributing graph is shown in 
Fig. |1|. Note that there are graphs that contribute to the diffractive gluon distribution 
at one lower order of perturbation theory, but they yield fewer powers of 1/xp. Thus we 
consider af contributions. 

We write the hadron momenta as 

p'X = (ri,o,o) 

p'a' = ((l-^-H'^(r^'^'-q) • (2.9) 

Thus the momentum transfer is|^ 

q'" = p'"A-PA' = ixpp\,-— ^^,q) ■ (2.10) 

/(^i — Xp)Pa 



' Note that in our notation, q^^ is not the virtual photon momentum associated with deeply inelastic 



scattering in Eq. (1.1), as is common in the literature. This should cause no confusion since in this section 



there is no virtual photon. 



The invariant momentum transfer is t = q^q^ 
the Xp — > hmit, we use the approximation t 



— q^/(l — Xp). Since we are interested in 
— q^. We suppose that q^ < M^. 



P ^pP. 




Figure 1: A particular graph that contributes to the diffractive gluon distribution. 

It is useful to work in a frame in which p^ ~ M, so that the initial hadron is approxi- 
mately at rest. 

We have a final state particle with momentum 

k2 



k^" 



■2k- 



k-,k) . 



(2.11) 



Since plus momentum g"*" = Xpp\ is delivered from the scattered hadron and plus momen- 
tum f3xjpp\ is removed by the measurement operator, plus momentum 



2k- 



[1 - f3)xpp\ 



(2.12) 



remains for the final state particle. Let us assume (as we will find, self-consistently, in 
this calculation) that there are no important integration regions with k^ -C M^ or with 
k^ ^ M^. That is, k^ ~ M^ in the integration regions that give leading contributions. 
Then k~ must be large, k~ > M/xp. The observation that the final state parton has large 
minus momentum is crucial to the calculation. 

Our analysis is simplified if we choose a physical gauge. Since the gluon in the final 
state has large minus momentum, it is natural |T^ to use the null-plane gauge A' = 0. 
(We make a few remarks on the calculation in Feynman gauge in Appendix C) 

We have seen that the final state hadron has a minus component of momentum of order 
M while the final state parton has minus momentum of order M/xp ^ M. For the virtual 
particles, we divide the integration over minus momenta into regions i~ r^ M and i~ ^ M. 
In a general Feynman graph, it is far from easy to make this division, but at order a^ the 
situation turns out to be quite simple. 
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In particular, the heavy quarks in our model hadron must have £~ ~ M. In order to 
leave the hadron in a color singlet state, two gluons (at least) must attach to the heavy 
quarks. In order that the intermediate heavy quark lines not be far off shell, the minus 
momentum delivered by each of these gluons must not be large. Finally, the only sink 
for the large minus momentum carried by the final state gluon is the vertex representing 
the measurement operator F^^ ^ which can absorb large i~ since it is evaluated at a fixed 
value of plus position, x'^ = 0. With a little thought, one realizes that all of the remaining 
internal propagator lines must carry large minus momentum. 

We are thus led to the picture shown in Fig. |^. In the lower subgraph, all loop momenta 
have i~ ~ M. In the upper subgraph, all loop momenta have i~ ^ M. Two gluon lines 
with £~ ~ M communicate between the two subgraphs. An example of a graph that 
contributes to Fig. ^ is shown in Fig. |I[ 




Figure 2: Structure of the diffractive gluon distribution. In the amplitude, two gluons are 
exchanged, one gluon is absorbed by the measurement operator, and one gluon is emitted 
into the final state. The subgraphs u and L are evaluated at lowest order of perturbation 
theory. 



The definition ( |2.1| ) together with the structure represented by Fig. ^ lead to the fol- 
lowing expression for the diffractive gluon distribution: 



df^mpx,,x,,ci',M) 



(J/djSD (I/O 



iE 



d^s 



4^/ 



d^S 



d'^kdk+ 



167r2 2n/3xppJ 2 ^, J (27r)4 J (27r)4 J {2Tr)^2k+ 
xTr {ulik, g, ^Of Uo{k, g, sYf,} {2ir)6{k+ - (1 - (3)x,p\) 
1 -iD{s')a^i -iD{q + s')py iD{s)^p iD{q + s)sa 



„/2 



le 



{q + s'y — ie s"^ + ie {q + s)^ + ie 



xLo{q,s',e,eXb Lo{q,s,e,e'] 



pa 
cd 



(2.13) 



Here the {27r)5{k~^ — (1 — f5)xpp\) results from performing the integration over y in 
Eq. (p.l|). We note immediately that the integration over the plus momentum of the 



final state gluon gives 



{2'n)5{k+ - {I - ^)x^p\) = — . (2.14) 



{27T)2k+ ' ' ' ' '^ ""^^ 2{l -(3)x^p\ 

The functions Uq and Lq are the amputated Green functions represented in Fig. ^ The 
subscripts distinguish these functions from simpler functions u and L that are defined 
below in terms of Uq and Lq and appear in the final formula ( p.20|) for the diffractive gluon 



distribution. The function uq carries two transverse vector indices that are not shown. 
One is the index j carried by the operator F^^ in Eq. (|2.1|) , the other is the transverse 
polarization index of the final state gluon. It also carries two color indices: the color index 
a of F^^ and the color index of the final state gluon. The notation 1i{uIuq\ denotes a 
summation over these polarization and color indices. The indices carried by Uq that are 
displayed are the color and polarization indices for the exchanged gluons. The counting 
factor 1/4 associated with the gluon exchange accounts for the two ways for attaching the 
labels s' and q + s' to the two gluons to the right of the final state cut and the two ways 
for attaching the labels s and g + s to the two gluons to the left of the final state cut. 
The notation iD{tj^i,/[(i'^ + ie] denotes the propagator of a gluon in null-plane gauge. The 
Green functions Lq depend on the transverse polarization vectors e and e' of the initial 
and final state photons, respectively. There is a sum over the two choices for each of these 
polarizations. 

We can now make a number of approximations that simplify Eq. ( |2.13| ). These approx- 
imations will be discussed again in the following section, but we outline them here in order 
to exhibit their effect on the overall structure of Eq. ( [2.13 ). 



First, consider the momenta s, g — s, s', and g + s' of the exchanged gluons. We have 
taken |q| < M. Let us assume that (as we will find, self-consistently, in this calculation) that 
|s| and |s'| are of order M in the integration region that gives the leading contributions. 
Now s^ must be of order M'^ /k~ , where k~ is the minus momentum of the final state 
gluon, in order that the gluons in the upper subgraph are not too far off shell. As we 
have just seen, k~ ~ M/xp, so s^ ~ XpM. On the other hand, s~ must be of order 
^'^ lv\ ~ M in order that the partons in the lower subgraph are not too far off shell. 
Thus s'^ s~ ~ XpM'^ <C s^. For this reason, the propagator for the exchanged gluon 
with momentum s is iD^p/ s^s^ ~ — iD-yp/s^. A similar power counting shows that only 
the transverse part of the momentum in each of the other exchanged gluon propagators 
contributes in the 1/xp -^ oo limit. 

Second, s^, s'^ and g^, being of order XpM, are negligible compared to the plus mo- 
menta in the subgraphs Lq. Thus we replace s"*", s'~^ and g"*" by in the subgraphs Lq. 
Then the integrations over s~^ and s'~^ can be associated with the upper subgraphs uq. Sim- 
ilarly, s~ , s'~ , and g~, being of order M, are negligible compared to the minus momenta 
in the subgraphs Uq, which are of order M/xp. Thus we replace s~ , s'~ , and q~ by in 
the subgraphs u. Then the integrations over s" and s'" can be associated with the lower 
subgraphs Lq. 

Third, the lower subgraphs in Fig. |^ are proportional to color matrices Tr[tatfe] = 0ab- 



Thus we can replace 



Lo(g,.,5,eOg.^^^c /°^^^''''f^-- . (2.15) 



Fourth, we will find that in the end only the Qf^u parts of the propagators of the ex- 
changed gluons count. Then the Lorentz index structure in Fig. |^ is Uq Lq^s. Since the 
partons in the upper subgraphs have very large minus momenta, this becomes approxi- 
mately Uq~Lq^. 



With these changes, Eq. (|2.13|) becomes 



df^f^{Pxp,Xp,q^,M) _ 1 1 1 ^ f £s f £s' f d^k 



E 



dx,.dt 647r2 47r/?(l -/3)x2(p+)2 2 f^ J (27r)2 J (27r)2 J (27r)2 

'ds'^ + _ „ \ / rds^ 



xTr j U^;^ul{k,q,s')^^ j ij^Uo{k,q,s) 



1111 

X— r 



2 c'2 („ I a'\2 



s^ (q + s)^ s'"^ (q + S' 

ds'- Lo{q,s',e,eXb^\ ( fds' Lo{q, 3,6,6')^ 



-i/i^ N-r Ui7 "'yji'" p-^«) 



Let us define simplified functions L and u by 

'ds~ Lo{q, s,e,e' 



!£ m-1 =P^AHc,,s,s,e') (2.17) 



and 



•ds+ 



/ds 
— Uo{k,q,s)-- = lCA9!uil3Xci,s) , (2.18) 

where I is a unit matrix in color space. We insert these definitions into Eq. ( |2.16|) and use 

TrI = N^ -1 . (2.19) 

This leaves a factor Tr{'u"''M} in which the trace is now over transverse polarization indices 
but not color. 

These changes give 

<;5(/5x„Xr,q2,M) _ 1 Cl{N^-l)gtl^ f d-'s f d^ f d'^k 



dxpdt 647r2 47r/?(l - /5)x2 2 fp J (27r)2 J (27r)2 J (27r) 

xTr{Mt(/?,k,q,s')M(/5,k,q,s)} 
1111 

X- 



o (q + s)2 s'2 (q + s')2 
xL{q,s',e,e')L{q,s,6,6') . (2.20) 



This formula gives the basic structure of the answer for the matrix element ( |2.1| ). We will 
now examine this in detail. 



3. DifFractive parton distributions and null-plane field 
theory 

In this section, we use the formulation of QCD quantized on planes of equal light cone 
coordinates to analyze the structure depicted in Fig. 0. In doing so, we will derive explicit 
expressions for the subgraphs u and L. This style of analysis is perhaps less familiar than the 
approach using covariant Feynman graphs in Feynman gauge, but it expresses the physics 
of the process in configuration space in a more transparent fashion. For those readers who 
prefer a standard covariant calculation, we present some of the essential steps in such a 
calculation in Appendix C We carry out the calculation of this section for the specific case 
of the diffractive gluon distribution. Then, in Sec. 3.4, we assemble the complete result for 
both the gluon and quark distributions. 



3.1 The upper subgraph 

Consider the function m(/5, k, q, s) that appears in Eq. ( p.20| ) and is represented by the 
upper subgraph in Fig. ^. The partons in u move with very large minus momentum through 
the gluon field that accompanies the heavy quark state that is approximately at rest. Our 
analysis is designed to draw the consequences of this, concentrating on the development of 
the states in space and time. 

External field to represent the exchanged gluons 

It is convenient to replace the gluons coming from the lower subgraph by an external 
color field A'^{x). We thus consider the matrix element 

M = J dy-e'/'-^P^^y- {k,s\F^{0,y~,Oy^\0)^ . (3.1) 



Here F~^^ is the measurement operator ( |2.6| ) for the gluon distribution function and k, s are 
the momentum and spin of the final state gluon. The matrix element is evaluated in the 
presence of an external color field A. We expand A4 in perturbation theory and extract 
the term A^o,2 proportional to two powers of QgA and zero additional powers of Qs- The 
coefficient in this term is the Green function Uq in Eq. (|2.13|) : 

Here A is the Fourier transform of A: 

A^M = J dx e^'-A^^{x) . (3.3) 

Since A represents the gluons exchanged with the lower state and these gluons are in a 
color singlet state (c/. Eq. ( p.l5| )), we will replace 

A;{-s) Aliq + s) -. 4. j;^ A;{-s) A'siq + s) . (3.4) 
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After making this replacement, the arguments that led to Eq. ( 2.16| ) lead us to anticipate 



that A^o,2 becomes much simpler in the 1/xp — > cxd limit, taking the form 

—— Uo{k; Xpp\, 0, q; s+, 0, s)"" . (3.5) 

Zn 



Using Eq. ( p.l8|) , we identify the function u that appears in the final formula (|2.20|) for the 



diffractive gluon distribution. Thus 

xl C A g!u{(3, k,q,s) . (3.6) 

where I is a unit matrix in color space. Our aim in this section is to calculate 7Vfo,2 in the 
1/xp —>■ oo limit, then to use Eq. (|3.6D to extract u{j3, k, q, s). 



In order to better illustrate the physical principles involved and to give some indication 
of how the present calculation would work at higher orders, it is useful to generalize the 
problem that we attack. Let us therefore consider a matrix element 

M= I dy-e'^^^^AV- {k,, si, . . . , fc^, s^|F,(0, y', 0)+^|0)^ . (3.7) 

Here ki,Si are the momenta and spins of one or possibly more final state partons, with 
k~ ~ M/xp. The matrix element is evaluated in full QCD in the presence of an external 
color field A considered at all orders of perturbation theory. As a;^. — > 0, the momenta 
k~ become large. On the other hand, the external field A stays fixed. That is to say, the 
minus momenta of the quantum particles become large while the minus momenta of the 
gluons in the field produced by the diffracted hadron stay of order M. 

The eikonal line 

The operator F~^^ in Eq. (|2.6|) contains the exponential of a line integral of the color vec- 
tor potential, which now includes both the quantum potential A and the external potential 

A, 

E{0,y-,O)=exp(-ig r dx~[A^{0,x~,0)+AtiO,x-,0)]tc) . (3.8) 

This eikonal line operator produces the same effect as if there were a special color octet 
particle, S, with a propagator 

and an interaction vertex with the color field 

- igt^nf" (3.10) 
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with n-e = e"*". We can build such a particle into the theory. Let the particle be created with 
an operator Q\{y~) and destroyed with an operator Qa{y~)- The commutation relation is 
[Qa{y~),Ql{y~)] = Sat- The action 



/ dx Ql{x~ 

J ~oo 



dx 



Sba - 9[At{0, X", 0) + ^+(0, X-, 0)]{Qha 



Qa{x-) (3.11) 



will produce the desired propagator and vertices. We need a notation for the states. We 
use QliO)\0) = \£,b) and 

Ql{0)\ki,si,...,kN,SN) = \ki,Si,...,kN,SN,£,b) . (3.12) 

Thus we can replace the operator F^^ by 

Q^{y-)aOiiO,y-,0) , (3.13) 

where, using Eq. (^.61), 

Oiiy) = -t/3x^p\ Aaiyy - d^A^iy)^ . (3.14) 

When making this replacement, we include the S particle in the final state and add the 
extra action ( |3.11| ) to the action for QCD in an external color field. 

In subsequent equations, we do not explicitly indicate the color index, b, for the special 
eikonal particle and we write Q'^{y^)O^{0, y^, 0) for Q'^{y~)aOl{0, y~, 0). Indeed, the color 
indices for all of the partons are also left implicit. 

Note that treating the eikonal factor as being produced by a quantum particle with 
special properties is more than just a technical trick. In the experimental determination of 
the gluon distribution, there is a short distance interaction that scatters a gluon constituent 
of the hadron and produces a system of jets with very large i~. (This applies for either 
the inclusive or the diffractive F2, for Xp ^ 1 or Xp -^ 1.) If the color field of the hadron is 
too soft to resolve the internal structure of the jet system, then the jet system looks like a 
color octet particle with infinite i~. We approximate 

2M-lil + .. (-•«)M2rn" + ••■)- ^ (-i9%W (3.15) 

and arrive at the interactions of the special eikonal particle. This idealization is incorpo- 
rated into the definition [10, 14 1 of the MS gluon distribution function. Deviations from the 



idealization are accounted for in the perturbative calculation of the hard scattering matrix 
elements for the physical process. 

Taking the high energy limit 

Our problem is now to find the limiting behavior of 

M = I dy-e'f'^^P^^y'{k^,s,,...,kN,SN,£\Q\y-)O\0,y-,0)\0)A (3.16) 
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when Xjxp — > oo and all of the k^ tend to infinity like M/xp. We analyze this problem 
using approximations suggested by the "aligned jet" picture of small x deeply inelastic 



scattering |13|. Similar problems have been addressed in many papers over the past few 



years (see, for instance, Refs. |T3, |]T8[ and |T^); the treatment in ||T8[ is especially close to 
that given below. Here we note that essentially the same problem was solved in Ref . |2^ , in 
which the authors addressed deeply inelastic scattering producing a ^Ji pair in an external 
f/(l) field in the high energy limit k^ — >■ oo. We simply adapt the treatment of [Q into 
the problem at hand. 



Our analysis of M. is based on null-plane-quantized field theory, as in Ref. |Q. (How- 
ever, we use the theory quantized on planes of equal x~ = (a;° — x^)/a/2 instead of planes 
of equal (x° + x^)/\/2 used in Ref. [|^. This is appropriate to a system with large minus 
momentum.) In this formulation of the theory, the role of the hamiltonian is played by 
P"*", which is the generator of translations in x~ . We refer to this operator as H. In the 
problem at hand, H is the generator of x~ translations in full QCD in the presence of the 
external color field A. To start, let us change to the interaction picture based on full QCD 
without the external field as the base hamiltonian Hq and the interaction with the external 
field as the perturbation V. In this picture, we write A4 as 

M = Jdy-e'^'^^P^^y'{ki,s,,...,kN,SN,S\ 

xU{oo,y-)[Q\y-)O^{0,y-,0)]jU{y-,-oo)\0) . (3.17) 

Here a subscript I on an operator denotes the operator in the interaction picture specified 
above. The evolution operator U is 



U{t2, n) = T exp (-i r dz-Vi{z-)\ . (3.18) 



Now, within the approximations used here, the interaction V does not produce soft parton 
pairs from the vacuum and thus V^|0) ~ 0. (We discuss the approximations and their 
validity at the end of this subsection.) Thus, we can replace U{y~ , — cxo) by 1: 

M = j dy-e'P^^^^Ay'{k,,s,,...,k^,SM^S\ 

xf/(oo,r)[Q^(y^)(^'(0,l/-,0)]/|0) . (3.19) 

Since Vi{z~') vanishes when ^ = 0, it is substantially non-zero only for z'^ ~ 1/M. 
(That is, the exchanged gluons carry plus momenta of order M.) On the other hand, the 
integral in ( p.l9| ) extends over a much larger range, \y~\ ~ l/{xpM). Thus most of the 
contribution to M. comes from the regions ?/~ ^ 1/M and y^ <^ —1/M. In the region 
y~ <ti —1/M, we can approximate U{oo,y') by f/(oo, — oo). In the region y^ ^ 1/^, we 
can approximate U{oo,y~) by 1. Then 

M ^ f dy-e'^^^PAy~{k,,s,,...,kN,s^,£\ 

J —oo 

xU{oo,-oo)[Q\y-)O^{0,y-,0)]r\0) 
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Jo 
x[Q\y-)O^{0,y-,0)],\0) , (3.20) 

where we are allowed to set the integration endpoints to zero instead of, say, ±1/M because 
the difference is of order Xp compared to the integral. Now, adding and subtracting a term 
in the integral over y^ < 0, we obtain 

J —oo 

X [f/(oo, -oo) - l][gt(^-)O^(0, y-, 0)],|0) 

/oo I _ 

-oo 

x[gt(^-)O:'(0,y-,0)],|0) . (3.21) 

The second term here is proportional to 

5{f5xpp\ + Y.kt) , (3.22) 

which vanishes because all of the terms in the argument of the delta function are positive. 
Thus 

J —oo 

x[t/(oo,-oo)-i][gt(r)(^'(o,y-,o)],|o) X (i + o(xr)) . (3.23) 

We can understand Eq. ( p.23| ) as follows. First, the operator Q'^O creates a gluon and 
one of the special eikonal particles. Then this state evolves according to QCD, possibly 
evolving into a system with more partons. Since it has very large momentum in the minus 
direction, its evolution in y~ is slow (except for the inevitable ultraviolet renormalizations). 
At y" ~ 0, this system of quarks and gluons passes through the external field. After that, 
it continues its slow evolution. 



With a straightforward derivation (which is given in Ref. |5^ in the case of abelian 
gauge theory), one finds that in the high energy limit the interaction with the external field 
becomes a very simple operator, which, following ||2^, we denote by F: 



f/(oo,-cx)) ^F . (3.24) 

The action of F is simply to produce an eikonal phase for each parton while leaving its 
minus momentum and its transverse position unchanged. If the parton is at transverse 
position b when it passes through the external field, then the phase is 

F(b) =P exp|-z^ /" °°ciz-^+(0,z-,b)t„| . (3.25) 
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Here the color matrices t^ are the generators of SU(3) in the representation appropriate to 
the color of the parton and the V indicates path ordering of the color matrices. In the case 
of the special eikonal particle, b = 0. Then Eq. ( |3.23|) becomes 



M = f dy-e'^^^P^^y~{h,si,...,kM,SN,S\ 

J —oo 

x[F-l][Q\y-)O^i0,y-,0)]i\0)x{l + Oix,)) . (3.26) 



The approximation ( |3.24| ) is, in essence, very simple. For a scalar parton with large 



minus momentum k that absorbs a soft gluon with momentum q, the approximation is 

-zg 2A;%) ^ ^,„ ' . {-ig 2k- e+) = — ^ {-zg e+) . (3.27) 



{k + g)2 + ie 2k~q+ + ie g+ + ie 

For partons with spin 1/2 and 1, the diagrammatic derivation is similar, but with a little 
work required to deal with the numerator structure. There is, however, a substantial 
difficulty that is related to the question of which gluons are soft and which are large k~ 
partons. For us, the external field represents the soft gluons and the partons all have 
large k~ . However, really a soft external gluon can produce multiple soft quantum gluons 
with an interaction that does not have the eikonal form. The soft quantum gluons can 
interact somewhere else with the large k~ partons (with interactions that do have the 
eikonal form). Thus the present derivation works only when such interactions among soft 
gluons are neglected. The present derivation also neglects interactions with large transverse 
momentum quanta that affect ultraviolet renormalization. Fortunately, in our application 
we extract a result at the lowest nontrivial order of perturbation theory, where none of 
these complications arise. 

The operator F 

Let us specify in more detail the matrix elements of the operator F between parton 
states. First of all, there are separate factors for each parton: 

(/Cl,Si,. . .,kN,SN,S\F\pi,s[,. . .,PN,s']y,S) = 

{ki, si\F\pi, s[) ■■■ {kj^,SN\F\p]\f, s'j^) {£\F\£) + -peTmnta.tio'n.s , (3.28) 

where the notation " + permutations" indicates that we should match identical partons 
in all possible ways. For a single parton state specified at time a;~ = by its momentum 
k = {k~,k) and null-plane helicity s, we have 

{k,s\F\p, s') = Sss' (27r) 2k-S{k- - j9")F(p - k) , (3.29) 

where F{p — k) is the Fourier transform of F{p) defined in Eq. ( |3.25| ): 

F(p -k)= fdh e'(P-'^)-^F(b) . (3.30) 
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For the special particle £ we have 

(^|F|^) = F(0) , (3.31) 

thus giving us back the eikonal phase factor that was part of the definition of the measure- 
ment operator, with the lower limit on the y~ integration approximated by — cxd. 

The Born approximation 

We now revert to the lowest order of perturbation theory at which a leading 1/xp — >■ cxd 
contribution is obtained. We take the order g"^ contribution to F — 1. For the evolution 
of the partonic state between time y~ and time zero, and then from time zero to the final 
state at x" ^ oo we take order zero of QCD perturbation theory. Then there is but one 
gluon in the final state and we have 

xexp(?pV)(p,s',£:|gt(O)0^(O,O,O)|O) , (3.32) 

where the subscript g^ reminds us that we are to use the g^ term in the expansion of [F — 1] 
and where p^ is the free particle plus momentum, 

P*-^- (3.33) 



For the matrix element of [F — 1] we have, using Eqs. ( p.29[ ) and ( |3.31| ) 



(fc, s,8\[Y- 1],2 |p, s\ S) = 6ss' (277) 2p-5ik- - p-) [F(k - p) F(0)] ^^ . (3.34) 

Also, we perform the y~ integration to produce an energy denominator. Then 

^"•^-IW? [^"'-P)^Wl- fepu"pV(2^-) <'='-P--l^'''°-°-°)l°) ■ (3-35) 
It will be useful at this point to adopt the notation 

IJXpPJ + p7(2A; ) 

Here ip represents the wave function of the gluon state just before it interacts with the 
external field. In ip, p is the transverse momentum of the gluon in the intermediate state 
and i is its transverse polarization {i = (1,2)). Since the minus momentum of this gluon 
has been set to the minus momentum of the final state gluon, we have 

k^ 

k- = ^ . (3.37) 

2(l-/3)x,p+ ^ ^ 
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Thus 



M 



0,2 



(27r)2 



F(k-p)F(0) ^^P'%Kp) 



(3.38) 



Evaluation of ip 

With the replacement in Eq. ( p.37|) , Eq. (|3.36|) reads 

,.(,.p)^.jg, (r-";f'"'ir (3.39) 

XppJ /Sk^ + (1 - /5)p2 
To complete the evaluation of ^, we need to evaluate the operator matrix element 

(fc-,p,s,6|Oi(0)|0) , (3.40) 



where, here, we have made the color indices a, 6 explicit. Using Eq. ( p.l4 ), 0^{y) 
-i(3xpp\ Aa{yy - d^AaivY, we obtain 



(A;-, p, s, 6|O^(0)|0) = -1 6ab [(3xrp\ e^ik', p, s) + p^e+{k-, p, s] 



(3.41) 

The polarization vectors, for transverse polarization s = 1,2, are those appropriate to 
A~ = EfauEfe: 



Thus 



e^{k-,p,s) = 6^' 



{k-,p,s\O^m0) 



+ n ^ ^ s 2(l-^)Srpi , 

e+(A;-,p,.) = ^f= ^, "^^ f 



■ ( o + .7S 2(l-/5)x^J9i . ,^ 

-I f3xjpp\ 6^' + — ^ / ^ p^p' 



(3.42) 



k2 



— Z 



. XpPj^ f o-i 2 



k2 



(/^k^^^'^ + 2(1 - (3)p'p' 



Inserting this result into Eq. (|3.39|) , we have 

f3k'^6^' + 2(1 - /3)pV 



V''^(k,p) 



/3k2 + (1 - /3)p2 



(3.43) 



(3.44) 



Expansion in powers of the external field 



We now expand the eikonal phase factor F(k — p) -F(O) in Eq. p.38| ), picking out 
the contribution proportional to two powers of the external field. We choose to call the 
momentum labels of the fields q + s and —s, and we symmetrize over which label belongs 
to which field. This gives 



M 



g^ f d^q f d\s 



0,2 



2 J (27r)4 J {271) 



At{q + s)Ai{-s) 
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x\^pjs{'k,'k) 



itatb 



itbta 



s+ + ze — s+ + te 



2n6{q^ 



+^j,(k,k-q) 



itbta 



+ 



itatb 



2Tx5{q^ 



_s^ + It —s+ + le 
-V'js(k, k + s) tbta 27r(5(g+ + s+)27r(5(s+) 

-V'j.lk, k - q - s) tatb 2n5{q+ + s+)27r5(s+)| 



(3.45) 



Note the denominators l/[±s+ + ze], which arise from the path ordering instruction in 
Eq. (|3.25|). For instance, 



Jdx jdy 9{y' 



> X-) 6-^(9++*+)^ e''^y~ 



27i6{q^ 



s+ + ie 



(3.46) 



Color singlet simplification 



We can simphfy the resuh in Eq. ( p.45|) if we recall that in our problem we are to replace 
the classical field product A'^{q + s)Aa{—s) by the corresponding matrix element of the 
quantum fields. In this matrix element, the color field is in a color singlet configuration, 
since it results from the scattering of a meson that starts in a color singlet state and ends 
in a color singlet state. Thus we are entitled to make the replacement (|3.4| ), 



At A 



b ""^a 



Jba 



t Zl+ /! + 



Then we can evaluate 



tah ^ba 



C, 



N, 



c 



(3.47) 



(3.48) 



since the generator matrices ta are in the adjoint representation of SU(3). This provides a 
great simplification because 



+ 



s+ + le — s+ + le 



27r(5(s^ 



(3.49) 



Thus 



M 



0,2 



9l Ca 
2Nl-\ 



dq c? q fds d s 



(27r)3 J (27r)3 
x|?/'j,(k,k)+^j,(k,k-q) 

-i/jjs{k,k + s) -^j^(k,k-q-s) I 



^+(0, q- + s", q + s) ^+(0, -s", -s) 



(3.50) 



This result has the form anticipated in Eq. ( |3.6| ). We are thus able to extract the 
function u{P, k, q, s). We have 



u{p, k, q, s) = ^(k, k) + ^(k, k - q) - ^(k, k + s) - ^^(k, k - q - s) 



(3.51) 



Here ip is given in Eq. (|3.44| ). Our evaluation of m(/?, k, q, s) is thus complete. 



3.2 The exchanged gluons 

Each of the exchanged gluons appears between the upper subgraph m, in which the 
partons have large momentum in the minus direction, and the subgraph L, in which the 
momentum components in the plus direction and in the minus direction are of the same 
order. The factor for one of the gluons is, in an obvious notation, 

where we have denoted the gluon momentum, either q^+s^ or — s^, by l^. Here v^ = (1, 0, 0) 
so that V ■ I = l^ . The vector v^ appears because we are using A~ = gauge. 

In going to the high energy limit, we have already made some replacements that simplify 
Eq. ( p. 521 ). First, according to Eq. ( |3.50D , only the plus component of the external field 



appears. Thus the factor ( p.52|) becomes 



^ 2FF3P -^-+ i.v r • ^'-''^ 



Second, in Eq. ( |3.5CI| ), our exchanged gluon fields are evaluated with /+ = 0. Thus Eq. ( |3.52| ) 
becomes 

We can make one more simplification. Because of gauge invariance, when we sum over all 
ways of attaching our gluon to the quarks in the lower subgraph, we can drop the term 
proportional to l^C^. (Here we use the fact that the two gluons must be in a net color 
singlet state, so that we effectively have abelian Ward identities.) Thus, our factor ( |3.52| ) 
finally becomes 

U-^C+ . (3.55) 



This result is built into Eq. ( ^1201) 



3.3 The lower subgraph 

We now turn to the lower subgraph. Consider the matrix element 

M = [^ Us fd^x e^PA'-y-im- {0\TO'{y)O{x)\0)j^ . (3.56) 

The operator O creates the initial heavy quark-antiquark state and then the operator O' 
destroys the final heavy quark state: 

0{x) = eeg e^^(x)7'^^(x) 

0'{y) = eeQe'^^{y)Y^{y) . (3.57) 
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Here e^ and e' are the polarization vectors for the initial and final state photons and e cq 
is the coupling of the heavy quark to the photon. The initial and final momenta of the 
photon are 

p^A = (p1,o,o), 

Pa' = (pl'.TT^.-q) • (3-58) 

We expand A1 in perturbation theory and extract the term A^o.2 proportional to two 
powers of QgA and zero additional powers of Qg. The coefficient of this term is the Green 
function Lq of Eq. ( p.l3[ ) : 



1 fdpt f d'^s f d'^w ~, , „ ~, .. 
^0,2 = t; ^ 77TZU TKZuMsy.Aiwy^ 



2J2ttJ (27r)4 7 (27r)4 ' 'f 

x{2nyS\pA + s + w-pA') Loiq,s,e,€'Ycd ■ (3.59) 

When we use the color singlet nature of Lq and make the simplifications that result in the 
1/xp — >■ cx) limit, Eq. ( p^.l6D leads us to anticipate A^o,2 takes the form 

/ds~ 1 
^ ^^^^3Y^o(q;0,3-,s,e,e')eV • (3-60) 



Using Eq. ( |2.17|) , we identify the function L that appears in the formula ( p.20|) for the 



diffractive gluon distribution. Thus 

Mo,2 =2jj^2j^J^ ^(^^ 0' «)c" Mw\ 0, -q - s); p\L{q; s, e, e') . (3.61) 

Our aim in this section is to calculate A^o,2 in the 1/xp — > cxd limit and then to use 
Eq. ( |3.61| ) to extract L(q;s,e,e'). 



The matrix element A4 is evaluated in the presence of an external gluon field A, which 
we think of as created by the scattering in the upper subgraph that we analyzed earlier. 
In this subsection, it is convenient to use a reference frame in which the minus momenta of 
the partons in the upper subgraph, k~ ~ li^ / (xpp^) , are fixed to be of order M as Xj= ^ 0, 
so that the external field can be regarded as remaining fixed in the Xp —>■ limit. Then p^ 
is large: 

p^ ~ M/xp . (3.62) 

We are now ready for the evaluation of A4. We use the same methods that we used 
for the upper subgraph. As we found for the upper subgraph, the interaction with the 
external field can be approximated by the eikonal operator F that supplies for a parton at 
transverse position b a phase 

F{h)=Vexp^~zgJ^^dy+A-{y+,OMta] • (3.63) 
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Note that here the indices + and — are interchanged compared to what they were in 
our evaluation of the upper subgraph. We do not, however, make this change exphcit 
by introducing a new name for the operator F, the phase function F[b), or its Fourier 
transform F{q), 

F(q)= fdhe'^'-^Fih) . (3.64) 

In A^o,2, we have the order g"^ term in the perturbative expansion of F, with everything 
else evaluated at order zero in QCD interactions: 

dpA' 



M 



0,2 



XI 



(27r)- JO 
\0'{y)[¥],.0{xm 



dy+ dy- dy / dx+ dx' d^ e'^'^''^-'''^-'' 



(3.65) 



Equation ( |3.65| ) is actually the main result. The rest of the derivation amounts to some 
straightforward manipulations. We first insert intermediate states into the expression: 



M 



0,2 



dPA' 



x{2n 



2tt jo 
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dy+ I dy 
dzi 



2zi 



dvi 



dy 

dZ2 

2^ 



dx 



dr2 



dx 

dz[ 

2^ 



dx e*^-*''^"*^-^'^ 



dr[ 



dz'2 
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dr' 



x(0|O'(y)|2;ip+,r;,s;;4p+,r;,,4) 

X {z'iPa'^ r'l. -s'l; 4pam 4. 41 [Fjgz kiPA. ri. Si; Z2Pa^ ra, S2) 
X {ziPa^ ri, si; Z2p\, rs, S2\0{x)\Q) . 



(3.66) 



Here the state created by the operator O consists of an antiquark with momentum {zip\, ri) 
and spin si and a quark with momentum {z2p\, Y2) and spin S2- After the operator F acts, 
we have a similar state with the momenta and spins denoted with primes. Next, we perform 
the space integrals. The integrations over y^ and x"*" give k~ denominators. The other 
space integrals give delta functions that can be used to eliminate some of the momentum 
integrations: 



M 



0,2 



4(2 



I-kV J ^^ J z 1 - z) J J z'{\ - z') J 



X 



q 



r'^ + M^ (q + r')^ + M 



n -1 r 



+ ie 



r2 + M2 



M2 



n -1 



l-z) 



+ ie 



z' (1 - z') 

X {0\O'mz'p\„ r', s[; (1 - z')p\„ -q - r', s',) 

X {z'Pa'^ r', s[; (1 - z')p\,, -q - r', S2I [F]g2 |zp^, r, Si, (1 - z)p\, -r, S2). 
x{zp+,r,s,;il-z)pi,-v,S2\OmO) . 



(3.67) 



We introduce here a wave function $ defined in Appendix B, so that we can make the 
replacements 

(^p+,r,Si;(l-z)p+,-r,S2|0(O)|O) _ 

-[r2 + M^]/z - [r2 + M2]/(l - z) 
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e CQ ^z{l - z) ^z, r, M, e),,s, , (3.68) 

and 

{Q\0'mz'p\,, r^ s',- (1 - z')p\,, -q - r^ g^) 
q2 _ [r'2 + M2]/z' - [(q + r')^ + M'^]/{1 - z') 

e eg ^z'H - z') $+(/, r' + ^'q, M, e^^s', ■ (3.69) 

For the matrix element of F we can write 

{z'p\,, r', s\] (1 - z')p\,, -q - r', s'^\[¥]g2\zp\, r, Si; (1 - z)p\, -r, Sa) = 

5s'^s,5s'^s, {27,)2zp\5{z'p\, - zp\) (27r)2(l - z)p\5{{l - z')p\, - (1 - z)p\) 
x[F(r'-r)F(r-q-r')]32 . (3.70) 

With these substitutions we have 

x^\z, r' + ^q, M, £')s2si [i"(r' - r) F(r - q - v%2^{z, r, M, e),,,, . (3.71) 

In this equation, we expand F{y' — r) F(r — q — r') in powers of gA and keep the order g"^ 
terms. There are four contributions. Using Eq. (|3.61| ), we extract L: 



L(q,s,M,£,£') = ^^— ^ / — — / dz (3.72) 

4 7r J [Z TTJ^ JO 

xTr{[-<l>"f(z,r + s + zq,M,£') -^^{z,r-s- {1- z)q,M,e') 
+ ^\z,r + zq,M,e') + ^\z,r- (1 - ^) q, M,^')] Hz,r,M,e)] . 
We find in Appendix B that 

^(z,k,M,e) = ^- ^ \(1- z)e-ak-a- zk-ae-cr + iMe -o-] . (3.73) 

This completes our evaluation of -L(q, s, M, e, e'). 



3.4 The gluon and quark distributions 

We now recap the result for the diffractive gluon distribution and give its extension 
to the case of the diffractive quark distribution. Let us introduce a parton index a, with 
a = g,q. Let us define functions Ua in terms of the functions u of Eq. ( |2.18| ) as follows: 



Ua{x,, (3, q, s, s') = 4^^^!.^^ 2 / 7^^ Tr {ui{(3, k, q, s') «,(/5, k, q, s)) , (3.74) 



47r/5(l-/3)x2 J {2n 
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with Ua given by the hnear combination (|3.51 ) of wave functions 



UaiP, k, q, s) = ipaiK k) - ipaiK k + s) + ipaiK k - q) - ^„(k, k - q - s) . (3.75) 

The difference between the gluon and the quark cases is in the expressions for the color fac- 
tors Ca and the wave functions ipa- The color factor for gluons may be read from Eq. ( |2.2CI| ), 

Cs = CliN',-l) . (3.76) 

A similar calculation yields the color factor for quarks: 

c, = C^iV, . (3.77) 

The wave function for gluons is given in Eq. ( |3.44|) , 



, _ /3k^^-^- + 2(l-/3)p-p^- 



We find in Appendix A that the wave function for quarks is 



//3(l-/3)k2 

Then we may rewrite the overall structure (|2.2CI|) of the result in the following form, for 
both the gluon and the quark distributions: 

<;g(/?x^,x^,q^M) _ 1 1^^ f d^s 1 r <fs' 1 

dxrdt QAtt^2^^J (2 7r)2 s2(q + s)2 y (2 7r)2 s'2(q + s')^ 

xL{ci,s,M,e,e')Ua{xp,P,q,s,s')L{ci,s',M,€,e') , (3.80) 



with Ua given in Eq. ( g^TJD and L in Eq. §J^. 



4. Behavior of the difFractive Green functions 



The result ( |3.80[ ) for the difFractive parton distributions is given in terms of two quan- 



tities: the functions U and the functions L. The latter contain the dependence on the 
specific diffracted system. The functions U, on the other hand, are universal Green func- 
tions. They control the process of diffractive deeply inelastic scattering for any small-size 
hadronic system. In this section we examine some of their properties and present results 
from the numerical integration of Eq. ( |3.80|) . 
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4.1 Ultraviolet and infrared finiteness 



Observe first, in Eq. (|3.80|) , the factors 1/s^, l/(q + s)^ from the propagators for the 



exchanged gluons (and analogous factors with s ^ s'). The poles at s = and s = — q are 
cancelled partly by U and partly by L. From Eqs. ( p.75| ) and ( p.72| ) we see that Ua oc |s|, 



L (X |s| as s — i> 0. Analogous behavior is observed for the other poles in s and s'. 

The Green functions Ua are constructed from the linear combinations of wave func- 
tions (|3.75| ) by integrating over the s-channel transverse momentum k. Note that each 



of the terms in Eq. ( |3.75| ) would give rise to an ultraviolet-divergent integration over k 
in Eq. (|3.74|) , but that the bad behavior cancels among the terms. This can be seen by 
expanding Eq. ( |3.75| ) for k^ — i> oo. Both the leading and next-to- leading terms in the 



expansion vanish. The first nonvanishing contribution to Ua is proportional to the second 
derivative of the wave function ip. This goes like 1/k^ at large k^. The net contribution to 
Eq. (|3.74|) from the large k^ region is therefore of the type (i|k|/|k|^. The physical reason 



for this cancellation lies with the partons being at the same transverse position in the limit 
k^ -^ cxD. Since the net color of the state is zero, the coupling of the state to gluons vanishes 
in this limit. Similar power counting shows that all of the integrations are convergent in 
the ultraviolet. 

Since the integrals are convergent both in the ultraviolet and in the infrared, we conclude 
that the transverse momenta dominating the integrals are all of order M. 



4.2 Asymptotic behaviors of Ua 

The functions Ua contain the dependence on the longitudinal momentum fractions Xp, 
p. The dependence on Xp is given simply by the overall factor l/Xp that defines the leading 
1/xp — > CX3 power (at which level we are working). The dependence on f3 is in contrast 
nontrivial. Detailed plots of the f3 dependence that we find by numerically integrating 
Eqs. ( |3.80|) and ( p.74| ) will be given in the next subsection. The limiting behaviors of Ug 



and t/g as /3 — >■ and /? — >■ 1, on the other hand, can be obtained analytically. Consider 
/3 ^ 0. We may take the limit inside the integral in Eq. ( |3.74| ) and evaluate the wave 
functions ijja in Eqs. (ITTSI) , (|379|) for /5 -^ 0. We get 



IT) i j / fQ 1^2 

V^y(k,p) = -P^[l + 0(/5)] , ^,(k,p) = ^^p-a[l + 0(/3)] , /5«1. (4.1) 



By using the expressions ( |4.1| ) in Eqs. ( p.74| ), (|3.75|) it is straightforward to check that the 



integral in d^k is convergent. From Eq. ( p.74| ) and Eq. (41) we conclude that 



Ug(xp-\ U,ocp\ /5 < 1 . (4.2) 

The behavior ( |4.2| ) can be understood on general grounds. The emission of a soft vector 
quantum has a 1/(3 spectrum, while the emission of a soft fermion does not. The behavior 
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( [4.2|) holds for any value of the transferred momentum q. In the case q = 0, in particular, 
the trace and the integral in Eq. ( |3.74| ) can be evaluated analytically in a simple way. 
The coefficients of the leading (3^0 terms take a fairly simple form as functions of the 
transverse momenta s,s': 



^l\2' 



[/,(x^,Aq = 0,s,sO = ^'^^yl~ ^^ -J2s.s' In f^^±4. 1 +s^ 



4:7!-'^ Xp 



+ In 
and 



s — s 



/\2' 



-,'2 



In 



(s + s' 

c'2 



r\2' 



P 1 
+ ln 



s — S' 



'\2 



In 



s + s 



/\2' 



(^S — S 
c/2 



l\2' 



[1 + 0(/3)] , /3 « 1 (4.3) 



^.(..,Aq^O,s,sO^|^2s.sMn(|±|:i;) 



[1 + Om , /3 « 1 . (4.4) 



(See Appendix D for calculational details and general expressions for the functions Ua-) 

Consider now /?—*>!. It is convenient to switch to a new integration variable v in 
Eq. ( pTJD by setting 



k = Av , A 
The functions f/„ are rewritten in the form 






U. 



gtca 



" 47r/52x2 J (2 



/72^^'("^"'^ 



(4.5) 



(4.6) 



The wave functions ipa of Eqs. (|3.78| ),( pT9D are rewritten as 



V^? 



v2 5'^ + 2 p' p^' 

v2 + p2 



V"^ + p"^ 



(4.7) 



The limit /3 — > 1 may be obtained by taking A — > inside the integral (|4.6| ). Consider 
the leading term in the expansion of the functions Ua (Eq. ( |3.75| )) about A = 0. By using 
the expressions ([4.7|) we get 



Un 



S^i - 



v2 S'^ + 2 s^ s-'' v2 6'^ + 2q' q> v^ 6'^ + 2 (q^ + s* ) (q^' + s^)' 



v2 + s2 



v^ + q2 



v^ + (q + s) 



[1 + 0(A)] 
(4.^ 



and 



Un 



V^ + S^ 



s ■ a 



v^ + q" 



qa 



iq + sj 



(q + s)-a [1 + 0(A)] . (4.9) 



We notice that the (3 —>■ 1 behavior of the functions Ua is different according to whether the 
transferred momentum q is zero or finite. Consider q = 0. Then observe that in this case 
the terms of order A° in Ug, Uq vanish. (In particular, for the gluon this term vanishes after 
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integrating over the angle of s, so that 2 s* s-^ -^ s^ 6^^ .) By expanding further in powers of 
A, for the quark we find 



" -wj-Z _|_ cZ I -.T-Z _|_ O'^ 



V^ + S^ V V 



4v-s 2s2 \ ^^^ „ 



(q = 0) . 



(4.10) 



For the gluon, the q = case has an additional cancellation at order A^, so that the first 
surviving term is proportional to A^: 



Mg -^ 2 A' 



2 V v^' 



6'' 



2v*v-'' Av-siv's^ + s'v^] 



v^ + s^ 



[4 (v ■ s)2 - v2 s^ - v^] (v2 6'^ + 2 s* sJ 



v^ + s' 



2^3 



By doing the P -^ 1 power counting from Eqs. 
behaviors 

t/, ex (1 - /?)2 , t/,cx(l-/3)S 1-/?«1 (q = 0). 



(v2 + S2)2 



+ 0{X') , (q = 0) . (4.11) 

(|I:^),(^:TU[),(P:TD we obtain the 

(4.12) 



In the q = case one can also determine the first nonzero f3 ^ 1 coefficients in a simple 
way by using the expressions ( 4.1U ),( i.ll| ) in Eq. (O) and performing the trace and the 
integral in d'^v. This gives 



f/g(x^,/?,q = 0,s2,s'2 
2s2s'2 



gtCliN!~l)_^^_^^, 



X 



An'^Xp 
2 {s^f - 5 (s^f s'2 - 5 s2 (s'2)' + 2 {s'^f 



Gs^'s 



2^/2 



+ 



+o{i-py 



(S2 - S'2)^ 
fs2)' - 5 (S^f S'2 + 10 (s'f (S'^f - 5s2 (s'2)' + (s'^)^ 



52-s'2)^ 



>"fe. 



1-/5 < 1 



(4.13) 



and 



Ug{xp,P,q = 0,s,s') 



^.^g^AT, 



TT^ X 



(1-/5) 



(S ■ S'f (S2 + S^2) 
(S2-S'2)2 



S2s^2 
g2 _ g/2 



/^2 1 



2 (s • s') 

fs2 - S'2)2 



1^772+^(1-/5)^ 



1-/3 < 1 . (4.14) 



(Again, we refer the reader to Appendix D for more details. 



In the case q 7^ 0, on the other hand, there are no such cancellations as in Eqs. ( [4.10| ),( ^rTT 
The leading A° terms in Eqs. (|47^) , ( ^^ do not vanish. The /? — >■ 1 power counting thus 
yields a constant behavior for the functions Ua- 



Ug,Ug^{l-(3)\ 1-/3«1 (q^O) 



(4.15) 
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It appears that near (3 = 1 the diffractive distributions will have a nontrivial q^-dependence 
and will be largest at nonzero q^. 



4.3 Results from Monte Carlo integration 

We are now in a position to determine numerical results for the diffractive parton 
distributions. To this end, we set up a Monte Carlo integration to evaluate the integrals in 
Eqs. (|3:8q),(F^,(pT2D. 



It is convenient to define the rescaled diffractive distributions ha 

haiP,qyM') 



X 



IM' dffli 



a^ ef) a't dxp dt 



(4.16) 



In Fig. ^ we plot results for the [3 dependence of these distributions at different values of 



lo"- 



10 ^ ^ 



10 ^ r 



10 ^ r 



n r 




solid: q=0 \ 

dotted: q/M = 0.2 'A 
dashed: q/M=l V^ 

dotdash: q/M = 2 



h. 




10 ^ 10" 10" 



Figure 3: The (3 dependence of the gluon (above) and quark (below) diffractive distribu- 
tions for different values of q^ ^ |t|. The rescaled distributions ha are defined in Eq. (|4.16|) . 

q. To emphasize the behavior in the regions of small j3 and large (3 we make a logarithmic 
plot in the variable (3 /{I — (3). This behavior reflects the properties of the functions [/„ 
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discussed in the previous subsection. In particular, for small {3 we have 

hg0^l3-^, hq(x(3^ {(3^0) . (4.17) 

For large f3, both the gluon and quark distributions evaluated at any finite q have a constant 
behavior: 

hg,h,^{l-(3f (/5-^l, q^O) . (4.18) 

The distributions at q = 0, on the other hand, vanish in the (3-^1 limit, because of 
the cancellations in the leading /3 — * 1 coefficients of the functions Ug, Uq observed in the 
previous subsection: 

/i,oc (1-/3)2, /i,oc (1-/3)1 (/3^1, q = 0). (4.19) 



^ 



^ 





' 




' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' . 




S 




. 




■ \ 




^ = 0.1 : 


.03 


- 


\ 


- 








\ 








\ 


.02 


- 




\_ 








s 








\ 








V 


.01 


~ 




^ - ~ 


















0020 


_ 


/■ 


1-|S-0.01 J 




■ / 




\ 








\ 


0015 


_/ 




N ~ 




'. 




\ 








X 


0010 


~ 




\ ~ 








\ 


0005 


~ 




~~'^-" 








, 1 , , , , 1 , , , , 1 , , , , 1 , , . . ■ 



2 3 

q/M 



Figure 4: The q dependence of the quark diffractive distributions for small or intermediate 
(5 (above) and for large /3 (below). 

The quantity that is perhaps most interesting for experiment at present is the diffractive 
distribution integrated over |t| ~ q^ from up to a value of the order of the squared heavy- 
quark mass, M^. Therefore in the following we will also consider the integrated gluon and 
fiavor-singlet quark distributions, defined as follows 



G=M 



M2 



rfq^ 



dxp at —" Jo 



M^ 



dq' ^^ 



dxpdt 



(4.20) 
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where M 



xl/{a- e 



2 rA 



OtZ 



and the sum in E runs over j = {u, u, d, d, s, s}. 



Fig. 1^ shows that the asymptotic large-/? behavior is reached for rather small values 
of 1 — /?, roughly of order 10~^. Furthermore, it shows that the asymptotic constants are 
numerically small compared to the values of the distributions at intermediate f3. Corre- 
spondingly, the diffractive distributions fall off as one approaches the small {1 — jS) region. 
The fall-off region may be most relevant phenomenologically, because likely it is to this 
range of j3 values (rather than to the asymptotic region) that current experiments on 
diffraction are most sensitive. 

The results in Fig. || indicate that the form of the q-dependence of the distributions 
is different at small or intermediate f3 and at large f3. We illustrate this in Fig. ^ for the 
quark distribution. Qualitatively, the behavior is the same for the gluon. While for small 
and intermediate (3 one has a monotonic fall-off with q, at large (3 the distribution peaks 
at a finite value of q. The decrease seen as q ^ when 1 — /5 -C 1 arises because of the 
cancellations in the functions Ug and Ug. 



10^ 



10" 



10" 



10" 




Figure 5: The gluon (hg) and quark (hg) diffractive distributions versus j3 (hnear scale). 

We observe from Fig. |^ that the gluon distribution is much larger than the quark 
distribution. The different order of magnitude at intermediate values of (3 (say, about 
P ^ 1/2) is roughly accounted for by the color factors in Eqs. ( p.76| ),( p?77D , Cg/cg = Tl 12. 
To facilitate the comparison, in Fig. ^ we display hg and hg in the same graph. Here we 
use a linear scale for /3. We see that the gluon remains large compared to the quark even 
at large values of /3. 
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Figure 6: The gluon (hg) and quark (hg) diffractive distributions for different incoming 
hadronic states. The dashed curves (photon case) are as in Fig. ^. The dotted curves 
(scalar meson case) and the dotdashed curves (gaussian case) have been rescaled by an 
overall normalization factor independent oi f3 (the same factor for gluons and quarks). 



The results described above have been obtained for the diffractive scattering of a model 
vector meson made of a photon that couples to heavy quarks. To what extent do they 
depend on the assumed incoming state? To address empirically this question, we first 
consider the case of a scalar meson that couples to scalar quarks with an interaction Xcpqq 
(A being a coupling with dimension of mass). A calculation analogous to the one described 
in Sec. 3.3 for the photon shows that in the scalar case the wave function (|3.73| ) gets replaced 
by 



$(^"^'^'')(z,k,M) = A 



zil 



k2 + M2 



(4.21) 



Then, we also consider a gaussian model in which the incoming bound state is described 
by the wave function 



^(.s'^^^^)(^z,k,M) = -^ exp (-kVM 



(4.22) 



We compute the resulting diffractive distributions for these cases. In Fig. ^ we plot the 
results. In this figure, the normalization for each of the two models ( ^4.21|) and ( |4.22|) has 
been fixed so that it includes coupling and charge factors as well as an overall arbitrary 
numerical constant. This arbitrary normalization factor is independent of f3 and is the same 
for quarks and gluons. Aside from this absolute normalization, we see that the results for 
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the diffractive distributions are remarkably stable against variation of the hadronic wave 
function. In particular, the shape in (3 of both the gluon and the quark distribution as 
well as the relative size between the gluon and the quark is qualitatively very similar for 
all initial states. We take this as an indication that such features of the diffractive parton 
distributions may hold with more generality. They do not so much depend on the specific 
form of the incoming hadron wave function, but rather they result from the color and x~ 
ordering constraints that give rise to the convolution formula (|3.80|) and the expressions 
( p.75| )-( ^^ ) for the Green functions [/„. 



Having given results for the diffractive parton distributions, it is of interest to ask what 
are the typical transverse sizes that dominate the answer. A measure of the transverse 
separations between the two outgoing particles in the upper subgraph (see Fig. 2) may be 
obtained by looking at the relative transverse momentum in the final state, that is, the 
momentum k in Eq. ( p.74|) . In Fig. ^we plot the distribution in k^ for the gluon operator 
matrix element at some fixed values of j3 with q = 0. We observe that there is quite a large 
spread in k^ values. In the range of j3 values considered, the typical value of |k| does not 
depend much on (3 and is of order M, about |k| ~ 3M. This plot provides a numerical 
illustration of the remark made below Eq. ( p.l2| ) about the dominant integration regions in 
the upper subgraph. A qualitatively similar behavior holds for the quark matrix element. 
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Figure 7: The distribution in the relative transverse momentum k for the gluon case at 
various fixed values of (5 and q = 0. 
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5. Ultraviolet behavior and renormalization group 

To the order in a^ at which we have worked so far, the matrix elements do not have 
ultraviolet divergences. Correspondingly, the calculation that we have discussed does not 
describe scaling violation. When additional gluons are emitted from the top subgraph 
in Fig. 1, on the other hand, ultraviolet divergences arise. The renormalization of these 
divergences leads to the dependence of the diffractive parton distributions on a renormal- 
ization scale /i. (The scale at which parton distributions are renormalized is often called the 
factorization scale in applications). Because the factorization theorem applies, this scale 
dependence is governed by the same renormalization group evolution equations as in the 



case of the inclusive parton distributions P, |10|, ^ 



The higher order, ultraviolet divergent graphs are suppressed compared to the graphs 
considered so far by a factor a^ log(/i^/M^). When log(/i^/M^) is large, these contributions 
are important, and thus evolution is important. On the other hand, when /i is of the same 
order as the heavy quark mass M, the higher order contributions are small corrections 
to the graphs considered so far. Thus one may interpret the result given in the previous 
sections as a result for the diffractive parton distributions at a fixed scale of order /i^ ^ M^. 
Then the diffractive parton distributions at higher values of fi"^ are given by solving the 
evolution equations with the results of Eq. (|3.8CI| ) as a boundary condition. 



We will explore the results of evolution in the next section. 



6. Evolution and diffractive deeply inelastic scattering 

Until now, we have considered diffractive parton distributions in a model hadron with 
a transverse size 1/M that is as small as one likes. Then lowest order perturbation theory 
is applicable. We have found that in this context we can evaluate the diffractive parton 
distributions as exactly as numerical integration permits. In Sec. 4, we have investigated 
some of the properties of the results. Among other things, we have found that the results 
are rather insensitive to the precise wave function of the small-size hadron; if we change 
wave functions we get diffractive parton distribution functions that are qualitatively similar. 

Now suppose that one had available a hadron of adjustable size R. For a very small size, 
R ^ 1/Aqcd, the diffractive parton distributions at any /x^ > l/R^ must be similar to those 
obtained by evolution beginning with Eq. ( p.80|) with M = 1/R. Let the size now increase. 



Longer and longer distances are now allowed to contribute to the diffraction process. What 
would the answer look like when R ^ 1/(300 MeV)? In a perturbation expansion, the 
result would be completely dominated by the soft region k± ~ 300 MeV. Then one possible 
scenario is that the diffractive parton distributions are radically different from those for a 
small hadron. A different, but conceptually related, scenario is that, even in the soft region, 
perturbation theory gives a description that is not too far off. In this case the diffractive 
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parton distributions would be similar to what one gets from the calculation presented 
above using evolution from the scale M k, 300 MeV to the multi-GeV scale relevant for 
experiments. On the other hand, as the size of the hadronic system increases, we may 
hypothesize that nonperturbative dynamics sets in that reduces the infrared sensitivity 
suggested by the perturbative power counting. As we go to larger and larger sizes, then, 
the distance scales that dominate the diffraction process, rather than continuing to grow, 
stay of the order of some intermediate, semihard scale 1/Msh- This suggests a conceptually 
different scenario for the diffractive parton distributions, in which the contribution from 
hard physics is enhanced with respect to the contribution from soft physics. Under this 
hypothesis, the diffractive parton distributions would be similar to what one gets from the 
calculation presented above using evolution from the scale M ^ Msh to the multi-GeV 
scale. 

Although this hypothesis does not have a firm theoretical justification at present, there 
are some indications in its favor. A set of indications comes from lattice QCD. Lattice 



investigations of glueballs ||2^ suggest that the correlation length for a color singlet pair 
of gluons is on the order of 1 GeV, not 300 MeV. Another set of indications comes from 
recent experimental measurements on the Xp dependence in diffraction. Recall that in our 
model calculation the leading dependence is x~^, corresponding to a pomeron intercept 
ap{0) = 1. The inclusion of higher order corrections of the type a^ In (l/xp) would, 
roughly speaking, modify this into a structure of the form ap{0) ~ 1 + const, x a^. The 
value measured in diffractive deeply inelastic scattering |^, Q is ap{0) ~ 1.15. The first, 
very general, observation is that this result is not inconsistent with evaluating a^ at a 
relatively short distance scale in the above formula for ^^^(O). More specifically, it has been 
stressed [|, ^, ^ that the experimental value of ap{0) — 1 given above differs by a factor of 
2 from the corresponding value measured in soft hadron-hadron cross sections. This may 
be taken to suggest that semihard physics dominates the diffractive parton distributions. 

In this section, we explore the hypothesis of a semihard scale by comparing its predic- 
tions to results from diffractive deeply inelastic scattering from protons at HERA. To carry 
out this study, we choose a value for M in Eq. (|3.80| ) and take the scale dependence of 
the diffractive parton distributions to be that given by the (two-loop) evolution equations 
( |1.2|) with the results ( p.80|) as a boundary condition at /x = M. We expect a semihard 



scale to be of the order of one or two GeV. In what follows we set M = 1.5 GeV. Note 
however that in this study the value of M is to be regarded as a free parameter to be ad- 
justed phenomenologically. If, for instance, the data prefer a value M -C 1 GeV, this would 
indicate that the diffractive parton distributions are completely dominated by soft physics 
but that nevertheless the perturbative result is not far off. If the data do not agree with 
the predictions for any M, then we would learn that soft physics dominates and completely 
transforms the answer. We will comment later on what happens when we vary the value 
of M. 

The study of evolution may be done for the distributions at |t| ~ q^ fixed or for the 
distributions integrated over \t\. Here we consider the integrated case. We determine the 
gluon distribution and the (flavor-singlet) quark distribution at the initial scale n = M 
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Figure 8: Evolution of the gluon (above) and singlet quark (below) diffractive distributions. 
The integrated distributions G and S are defined in Eq. ( [4.20| ). We assume the initial scale 
to be M = 1.5 GeV and we use evolution equations in next-to-leading order (two-loop). 



by integrating the result ( p.80|) according to the definition (|4.2CI|) . We then compute the 
evolution of these initial distributions up to fi = Q for different values of Q. In Fig. |^ we 
show the results. 

It is interesting to compare directly the distributions that we obtain from this calculation 
with the ordinary (inclusive) parton distributions. In Figs. ^and|lO|we report the gluon and 
up quark distributions at a certain scale (Q^ = 20 GeV^) for, respectively, the diffractive 
case and the inclusive case. The diffractive result (Fig. |^) is from our calculation. The 
inclusive result 
proton PI 



I^Fig. |10D is from the standard set CTEQ4M of parton distributions in a 
In the diffractive case we look at the j3 dependence and plot /5^ times the 
distributions. In the inclusive case we look at the x dependence and plot x^ times the 
distributions. 

Recall that the absolute normalization of the diffractive distributions has been obtained 
by dividing out numerical factors and coupling factors associated with the incoming state 
(see Eq. ( [4.20|) ). Therefore the absolute scale of Fig. ^ compared to Fig. |10| has to be 
regarded as arbitrary. In contrast, the relative normalization between gluon and quark as 
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well as the shape in /3 is determined by our calculation. Fig. 10 shows that, for ordinary 
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Figure 9: The P dependence of the gluon and up quark diffractive distributions (multiplied 
by /5^) at Q^ = 20GeV^. The overall normalization constant C is arbitrary. 
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Figure 10: The x dependence of the gluon and up quark inclusive distributions (multiplied 
by x^) at Q2 = 20 GeV^ (from the set CTEQ4M). 
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partons in a proton at the Q^ scale considered, the gluon is dominant for small momentum 
fractions, but, as the momentum fraction increases to about x ~ 0.2, the up quark starts to 
dominate over the gluon. The behavior changes dramatically in the diffractive case (Fig. ^. 
The gluon distribution is broader and stays large compared to the quark throughout the 
range of momentum fractions. The origin of this is in the behavior observed at a fixed scale 
M (see Figs. ||, H). Fig. | illustrates that the feature found in the fixed scale calculation 
persists qualitatively after the inclusion of loop corrections through perturbative evolution. 
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Figure 11: Scaling violation in the (flavor-singlet) quark distribution S at moderate values 
of momentum fractions. Above is the case of the diffractive distribution, below is the case 
of the inclusive distribution (from the set CTEQ4M). 



This behavior has consequences on the pattern of the scaling violation in Q^. As can 
be seen from Fig. ^ the diffractive distributions grow with Q"^ at low (3 and decrease with 
Q^ at high (3. In particular, for the quark the stability point at which the behavior changes 
is about /3 ~ 0.5. This is to be contrasted with the case of the ordinary (inclusive) quark 
distribution in a proton, for which the stability point is at x ~ 0.08. There is a large, 
important range of moderate values of momentum fractions, say, approximately, 0.1 to 0.5, 
in which the ordinary quark distribution is flat or weakly decreasing with Q^, while the 
diffractive distribution is rising with Q^. This is illustrated in Fig. |ll[ The explanation 
for the rise in the diffractive case lies with the gluon distribution being dominant even at 
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large momentum fractions (Figs. ^Q). As Q^ increases, gluons splitting into qq pairs feed 
the quark distribution and cause it to grow in the region of moderately large (3. 



Let us turn to the calculation of the deeply inelastic structure function F2 ■ Having 
determined the diffractive parton distributions and their evolution, we may use the fac- 
torization formula ( |1 . Ij) to compute results for F2^^. We evaluate -^2^'''^ in next-to-leading 



order by using the one-loop expressions [0 for the hard scattering functions Fa- Note that 
the use of the factorization theorem allows us to systematically take into account correc- 
tions to the diffractive scattering beyond the leading logarithms for both the quark and 
the gluon contributions, much as in the case of inclusive deeply inelastic scattering. We set 
the factorization scale /i^ in Eq. ( |1 . 1| ) equal to Q^. In Fig. ^ we plot the results for F2^^ 
versus P at different values of Q^ 
of Eq. (O) by 



Here F^^^ is related to the differential structure function 



p,m 



<Ljr 



9 4 



«„ 



M^ 



d\t\ 



dFf^ 

iJjJbp (JiL 



(6.1) 



where a is an arbitrary numerical normalization. In Fig. |T^ we show these results along 
with the ZEUS data for the proton diffractive structure function at the same values of 
Q^ [0 . These data are obtained from the structure function integrated over t by fitting the 



Jbj. 



dependence to a power x^^"'' and extracting the coefficient of x^^"^. 
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Figure 12: The [3 dependence of the diffractive structure function F2^^ for different values 
of Q^. We compute F2^^ in next-to-leading order. 



Notice the main qualitative features of the curves in Fig. 12. For 0.2 < j3 < 0.8 the 



diffractive structure function is rather flat in (3. As Q increases, F2 increases for (3 < 0.5 
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Figure 13: Same as in Fig. |T2|. Also shown are the ZEUS data from Ref. 



reflecting the behavior already observed for the quark distribution S. For 0.2 < f3 < 0.8 
the Q^ dependence is quite mild. Similar rather flat dependences on both /3 and Q^ 



are 



striking features of the data for -F^'^ and distinguish the diffractive -F^'^ sharply from the 
inclusive Fo. 



We observe from Fig. |T3| that these features are in qualitative agreement with what is 
seen in the HERA data on proton diffraction (except for the two data points at the smallest 
values of /3: see below). Recall that, given the value of M, both the Q^ dependence and the f3 
dependence of F2^^ are determined from theory. Only the overall normalization, associated 
with the absolute normalization of the diffractive parton distributions, is free. (The t 
dependence is also determined in principle. In the results presented here t is integrated 
over.) The agreement represented in Fig. |1^ between the predictions and the data is 
evidently not perfect, but given the simple nature of the theoretical calculation, one may 
suspect that this agreement is telling us something. 



Recall that the predictions represented in Fig. |I3| are based on the choice M = 1.5 GeV. 
If we decrease M to 1 GeV, we find that the agreement between the predictions and 
experiment is also not too bad, but if we choose M much below 1 GeV then we find that 
there is too strong a dependence on Q^ and too much of a slope in /?. Thus a semihard 
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scale for M seems to be preferred by the data. 



In the region of small /5, the curves of Figs. |T2|,|T3| have a different behavior from that 
suggested by the two data points at the lowest values of [3 and lowest values of Q^ {Q^ = 
8 GeV^ and Q^ = 14 GeV^). If further data were to confirm this difference, this could point 
to interesting effects. Here we limit ourselves to a few qualitative remarks. As far as the 
theoretical curves are concerned, we note that the diffractive distributions that serve as a 
starting point for the evolution are fairly mild a.s (3 —>■ 0. The gluon distribution goes like 
1//3, while the quark distribution goes like a constant (see Eq. ( |4.17| )). The small- /3 rise 
of the structure function ^2^'''^ in the curves of Figs. p!2| , p!3| is essentially due to the form of 
the perturbative evolution kernels. As regards the data, it has been observed |2^ that for 
small (3 the experimental identification of the rapidity gap signal may be complicated by 
the presence of low p± particles in the final state. If the current data hold up and especially 
if the same features are observed at lower values of (3, it would be interesting to see whether 



detailed models for the saturation of the unitarity bound |^ , which might bring in a new 



kind of physics for very small Xp, could accommodate this small (3 behavior. 



7. Conclusions 

We have analyzed diffractive deeply inelastic scattering in QCD, drawing on two main 
ideas. One is the notion of factorization of the hard scattering. For processes with only 
one hadron in the initial state, this allows us to introduce diffractive parton distributions, 
defined through certain measurement operators in terms of the fundamental quark and 
gluon fields. The other is the notion, widely used in small x physics, that the spacetime 
structure of the diffraction process looks simple in a reference frame in which the struck 
hadron is at rest (or, at least, has small momentum). That is, the process is dominated by 
configurations in which the parton probed by the measurement operator is created at light 
cone times x~ far in the past and much later, following a slow evolution in x~ , interacts 
with the color field of the incoming hadron. 

We find that this physics can be described most transparently by using a hamiltonian 
formulation in which the theory is quantized on planes of equal light cone coordinates. We 
develop this description in detail along the lines of Ref. l2( 



Using this method, we examine the problem of diffractive deeply inelastic scattering 
from a color source with small transverse size 1/M. Because the size is small, perturbation 
theory is applicable. Using perturbation theory at lowest order (order a^) for the diffractive 
parton distributions, we find a solution in closed form. This solution has the property that 
(except for the overall normalization factor) it is independent of the size of the source. It 
is also scale invariant, since the diagrams are ultraviolet convergent. The scale dependence 
that arises from diagrams with more loops is to be included by solving the renormalization 
group evolution equations with the results of the lowest order calculation as input at yU = M. 
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We find that, although the diffractive distributions are generally softened by evolution, 
certain distinctive qualitative features survive. Both the diffractive quark distribution and 
the diffractive gluon distributions fall off less quickly as (3 increases than the comparable 
inclusive distributions fall off as x increases. Furthermore, the diffractive gluon distribu- 
tion is much larger than the diffractive quark distributions. This has implications for the 
diffractive structure function: ^2^^^ is rather flat as a function of both (3 and Q^ as long as 
Q^ is not too much larger than M^. 



We observe that these qualitative features can be seen in the data from HERA on 
the proton diffractive structure function ^2^'''^ in the range 0.2 < (3 < 0.8 and 8 GeV^ < 
Q^ < 60 GeV^. If the small hadron calculation is regarded as a model with an adjustable 
parameter M, then we find rough agreement between the model and the experimental 
results if we take M to be Msh ~ 1-5 GeV. We conclude that the model may not be too 
far from reality provided that some essentially nonperturbative effect intervenes to provide 
a semihard scale Msh in diffractive deeply inelastic scattering. 
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Appendix A. Wave function for quarks 

In this appendix, we evaluate the wave function ipssi}^-, p) for the antiquark state created 
by the quark measurement operator. We follow the gluon case closely, merely noting what 
needs to be changed to deal with the quark measurement operator instead of the gluon 
measurement operator. 

In Eq. ( |2.7| ) we have the operator 

^^ ^(0) Et(o) . . . ^+ ... E(0, y-, 0) vI/(0, y-, 0) 
= ^^^\Q)E\Q)... i7-7+ ■■■E(0,r,0)vl/(0,r,0) . (A.l) 
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We use a representation of the 7 matrices with 

Then 

47V = (J °) . (A.3) 

Let us break the Dirac field into two parts: 



^ 



L 



Then the operator in Eq. (2.7) is 



^^ vl/t^(0) E^iO) . . . EiO, y-, 0) ^uiO, y', 0) . (A.5) 

Thus in order to compute the upper subgraph for quarks, we need to compute a matrix 
element analogous to that in Eq. ( p.l|) : 



M = J dy-e'^^-PAy'[px,py^]y^{k,s\E{0,y-,O)%{0,y-,O)\0)^ . (A.6) 

The Dirac field operator ip creates the antiquark state that, after interaction with the 
external field A, becomes the antiquark state {k,s\. Proceeding as in the gluon case, we 
obtain the analogue of Eq. ( p.32| ), 

M = J^J^Jd'pJ2{k-,^;s■,S\[F-l]^.\p-,p■,s'■,S)^|J,A^,p) , (A.7) 



where the wave function analogous to that in Eq. ( ^.361 ) is 



J, rw nl - , l^^^-^i V^' {k-,p;s'\%iO)\0) 



with k = k^/[2(l — P)xpp\] as in Eq. ( p.37|) . Here s, s' and s are two component spinor 



indices. We evaluate the matrix element in free field theory, so the result is simple: 

{k-,p,s'\%{0m = Vlj{k-,p,s') , (A.9) 

where Vu{k~, p, s) is a two-component spinor consisting of the top two components of the 
four-component Dirac spinor for antiquarks, V{k~ ,p,s): 

^(^-.-)-(v:,t':pM)- f^-^") 

It remains to specify precisely the spin, and here we come to an important technical 
point. As we have seen, the operator that we need is expressed very simply in terms of the 
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upper two components of the Dirac field. These are the components that have a simple 
partonic interpretation for a system of partons with large momentum in the plus direction, 
the direction in which the hadron is moving. However, our derivation has been based on 
null planes with fixed x~ , which is the natural formulation of the theory for a system of 
partons with large momentum in the minus direction. In this formulation of QCD, it is the 
lower two components of the Dirac field that are the independent degrees of freedom for 
the quarks. At any given x~ , the upper two components are given in terms of the lower two 



at the same x by an equation of constraint [[Iq , ^ . The lower part of the Dirac field has 



two components, which create the two antiquark states with the corresponding null-plane 
helicities: 

(fc-,p,S2m(0)|0)oc5,,,, . (A.ll) 

With this definition of helicity, the antiquark helicity is preserved as the fast antiquark 
passes through the external field. 

We conclude that with the appropriate definition of the antiquark spin the lower two 
components of the spinors V are simple: 

Vl{k-,p,s') = [V2k~Y/^S-s,, . (A. 12) 

Here the normalization is fixed to give the conventional normalization for the four- component 
spinors, V{k, s)'y^V{k, s) = 2k^. The upper two components of the spinor are related to 
the lower two components by the free Dirac equation 

Vt/(A;-,p,s) = -^p-aVL(A;-,p,s) . (A.13) 

Thus 

(fc-,p,g^|0^"(O)|O)= ~' p-a,v . (A.14) 

V2k-^ 

We now insert this result into Eq. (|3.39| ) in order to obtain the wave function for the 
antiquark state created by the measurement operator. For our application in Eq. ( |A.7| ), we 
want s' = s because of null-plane helicity conservation and we want k~ = k"^ / [2{1 — f3)xpp\]. 
With these replacements, we find 

, „ . 1,3(1 -/3)k2]'" ,. ,^, 



This expression for ipss is ready to be inserted into Eq. (|3.75|) for the upper subgraph U. 
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Appendix B. Wave function for the incoming heavy- 
quark system 



In this appendix, we evaluate the function $ used in Eqs. (|3.68|) , (|3.69|) . We define $ as 

{zp\, n, gi; (1 - z)p\, ra, 52^^^(0)7^^(0)10) _ 
(n + V2Y - [rl + M^]/z - [rl + M2]/(l - z) 



z{l-z)^{z,{l-z)Yi-zY2,M,e)s,s^ . (B.l) 

It is not immediately obvious that $ depends on the combination (1 — 2;)ri — ZT2 of the 
transverse momenta and not on ri and Y2 separately. This is a consequence of the invariance 
of the wave function under the subgroup of the Lorentz group in which ri — i> ri + z\ and 
i"2 ^ i"2 + (1 — z)w pO]. We shall find this property by explicit calculation. 



First, the denominator is 

The numerator is 

£{p\, ri + r2)^ W((l - z)p\, rs, S2) 1^ V{zp\, n, Si) . (B.3) 



For the Dirac algebra, we use the representation ( |A.2| ) of the Dirac matrices. Then, using 



e^ = Q gauge for the polarization vector, the numerator has the form 



t ;,tN ( V2^ ^^± ■^±\ f^u 



We need the quark and antiquark spinors lA and V. In our application, the quark and 
antiquark move through an external field with a large momentum in the plus direction. The 
external field leaves the spin unchanged provided that we use null-plane helicity adapted 
to the plus direction. Then the upper two components of U and V are simple: 



-'SS2 1 



U^{{l-z)p\,r2,S2) = [V2 {1 - z)p\f' 6s 

Vlj{zp\,T,,S2) = [V2zp\Y/'6ss, . (B.5) 

The lower two components are then given by the free Dirac equation: 

UI{{1-z)pU2,S2) = ^|-^^W^((l-;.)pi,r2,.2) , 

V2(l - z)p^ 

Vi{zp\,v„S2) = '^^-^-^Vtj{zp\,v,,S2) . (B.6) 

Finally, we need the polarization vector e{j>\, ri + Y2Y of the photon. Since p^e{j)Y = 
in £+ = gauge, we have 

e- = i^l±lAll . (B.7) 

Pa 
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With these ingredients, we obtain Eq. (p.l| ) with 

$(;2,k,M,e) = ^- ^ \(l- z)£-a'k-a- zk-ae-a + iMe-a] . (B.8) 



Appendix C. Covariant formulation 

In the main text we have derived the results for the diffractive parton distributions 
collected in Sec. 3.4 by using the null-plane formulation of perturbation theory and empha- 
sizing the picture of diffraction scattering in configuration space. Here we outline the main 
aspects of the alternative derivation based on the covariant formalism in momentum space. 
We limit ourselves to highlighting the most important ingredients of the calculation. 

We first discuss the gauge for the gluon field. The conceptually simplest choice would 
be to use A~ = gauge. However, some calculational simplicity can be achieved by using 
Feynman gauge. It is necessary only to replace the sum, —g'^", over four polarizations for 
the gluon that enters the final state by a sum 

J2e'^{k,X)e''{k,X) (C.l) 

A=l 

over physical polarizations, with the choice e~{k, A) = for the polarization vectors. Recall 
that 

~g>^^ = Y.e^{kA)e^{kA)- \^ , (C.2) 

where n'^ is the gauge fixing vector defined hy n ■ A = A^ . The unphysical polarizations, 
represented by the second term, do not contribute to the result after summing over graphs. 
Thus we can drop them. Effectively, then, we have A~ = gauge for the final state gluon. 
We can retain the simple Feynman gauge propagator, —ig^^/{k'^ + ie) for virtual gluons. 

Let us now see which graphs contribute to the result in the 1/xp — > oo limit. Consider 
the upper subgraph in Fig. 2. The definition of the Green function Ua contains integrations 
over the plus momenta s^ ,s'^ (see Eqs. ( p.l8|) , (|3.74|) ). We recognize that only certain 
topologies give rise to pinch singularities in the complex s^ plane. Consider for instance 
the two contributions in Fig. |14[ From the eikonal propagator and the t-channel propagator 
we get, for graph (a), 

t(") = / — - (C 'W 

+ 7-00 2t^ [s+ + ie\[{k-q-sf + ie\ ^'' 



and for graph (b) 



7-(^) = / i! (C A\ 

+ 7-00 2ti [s+ + k+ + ie\[{k + sf + ie\ ' ^ ' 



where N^"'^ and A^^'') denote spin numerators. Using the kinematic relation (|2.12 ) and 
the fact that in the hmit l/xp — »• cxd we only need to evaluate the upper subgraph at 
s~ = q~ = (see Sec. 2.2), the position of the poles in the two cases is as in Fig. |15[ 
In the case of graph (a) the poles pinch the contour and one gets a leading 1/xp -^ oo 
contribution, while in the case of graph (b) one may deform the contour away from the 
poles and neglect this graph as l/xp -^ oo. 
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Figure 14: Two examples of Feynman graphs contributing to the amplitude for the gluon 



Green function Ug. 



The gauge choice e~ = for the final state gluon is crucial for this argument to go 
through in the form described above. With this choice, the spin numerators behave like 
(s"^)*^. Then the integrals ( |C.3|) , (|C.4|) are sufficiently convergent at large s^ to allow the 
contour deformation. 
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Figure 15: Poles in the complex s^ plane from the graphs of Fig. 114 



Consider now the two graphs shown in Fig. |TB|, in which a gluon is emitted into the 
final state from the eikonal line. The contribution in Fig. p!6|(a) is disallowed by color 
conservation, because the t-channel gluons are in a color-singlet configuration while the 
eikonal line carries color octet charge. The contribution in Fig. |16](b) may be dealt with 
by s^ contour deformation, by an argument similar to that given above. 

Having seen why certain graphs do not contribute in the 1/xp -^ oo limit, we turn to 
the calculation of the result for a graph that does contribute. Consider the contribution 
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Figure 16: Gluon emission from the eikonal line. 



to the gluon Green function Ug from the product of the amplitude in Fig. 0(a) with the 
complex conjugate amplitude of the same form. This contribution can be written as 



2ti (3Xf 



27r 



27r 



(2 7r) 



X 



-F^"^ M^^ F"^ 



[s+ + ie] [{k-q- sy + i e] [{k - q - s'^ - i e] [s'+ - i e] 



.(C.5) 



where F'^^ is the gluon operator vertex, M^j, is the polarization tensor associated with the 
squared amplitude, and the resulting expression for the spin numerator is 



-F^^M^^F''^ = A{k-fU(3xrpf) - 



2 (3xpp\ 



k' 



(k-q-s)2 + (k-q 



[(k-q-s)-(k-q-sOf 
{k~Y 



(C.6) 



We may perform the integrations over the plus and minus components of the momentum 
k by using Eq. ( |2.14| ). The s"'',s'"'" integrations are of the form (|C.3| ). By carrying out this 
integral we get 

(C.7) 



I[ 



(a) 



1 



/3kV(l-/3) + (k-q-s)2 



Then we obtain 



C/f(x„/3,q,s,s') 



9tCliN^-l) f rf^k 



Cp 



J ^,{2mr+2(^{i-(3)\^ 



+ (k-q 



,/\2 



ATTp{l-p)xj 

4(l-/3)M(k-q-s)-(k-q-sO]'} 



q 



/{[/3k2 + (l-/3)(k-q-s)2] [/3k2 + (l-/5)(k-q-s'f]} 



(C.^ 



We recognize in this expression the contribution to Ug from the last term in Eq. (|3.75|) . 
Indeed, the expression (|C.8|) can be recast in the form 

U!^'\xp,p,ci,s,s') = llf^f_^-^^ ^ , 

(C.9) 
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with ip being the wave function given in Eq. ( p.78 ). 



The other terms in Eq. ( p.75|) are obtained in a similar way by adding the graphs in 



which the exchanged gluons attach to the same hne and using the color singlet projection 



Appendix D. An integral representation for the Green 
functions Ua 

In this appendix we give an integral representation for the Green functions Ua alternative 
to the one given in the text, Eq. ( |3.74j) . It is obtained from Eq. ( |3.74|) by performing the 



trace, introducing a Feynman parameterization to combine the denominators and carrying 
out the integral over k. The resulting formulas are lengthy and not as compact as the 
formulas used in the text. They have the advantages, though, that the spin structure is 
explicitly worked out and that the cancellation of the singular terms at large k^ is also 
explicitly worked out. They may be useful for numerical evaluation. 

The functions Ua are written as 



16 TT^ [3 Xp 

X [:Fa,j In Ml^ + %f-] , (D.l) 

a,ij ) 



M 



where the arguments (q, s, s , /3, x) of the auxiliary functions JF, ^, M. are suppressed. The 
color factors Ca are given in Eqs. ( p.76| ),( |3r77D . The summation labels i, j run over the four 
momentum configurations appearing in the definition of Mq (see Eq. (|3.75| )) 

u,(/3,k,q,s) = ^(-l)Va(k,k + r») , 

i=0 

where for i, j = 0, 1, 2, 3 r^ ,r2 have the values 

r? = {0,s, -q, -q-s} , r^-'^ = {0, s', -q, -q - s'} . 

The expressions for TaM-, Ga,ij and 7V1?- are 



or 



^a,ij = fa,i{x, 7)r? ■ rP + /a,2(a;, 7)rP ■ r^-'^ + (x -^ (1 - x), vf -^ r^ 

+^,,4(x,7)(rP ■ r?)(r« ■ r«) + (x ^ (1 -x),r« ^ r^^'^ 
Ml = 7x(l - 7x)rf^ ■ rP - -f^x{l - x)rP ■ 4^'^ + (x ^ (1 - x), r^ -^ 4^'^ 
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where 7 



(3. The new auxihary functions / and g are simple polynomials in x: 



/q.2 = 

9q.l = 

9q,2 = 

9q,2 = 

9qA = 

in the case of quarks and 

fg,i = 

fg,2 = 

9g,l = 

9g,2 = 



9g,3 



9gA 



4(1 - 7)a;(2 - 37a;) 

((1 -7)(-l+ 27- 1272(1 -x)x))/7 

-2(l-7)7V(l-7x) 

-(1 - 7)7^(1 - x)x[l - 27(1 - x)x] 

2(1 - 7)7(1 - x)x{l - 7 + 27^(1 - x)x) 

2(1 - 7)7x^(1 - 7(3 - 2x) + 47^(1 - x)x) 



(-2(1 - 2(1 + 27 + 3-f^)x + 6(7 + 7^)a;'))/7 
2(-l + 27 - 6(1 + -f^)x{l - x)) 
27x^(1 + 7(-l + x)2 - 27^ + 7^a;2) 
7([l - 2(1 + -f^)x{l -x) + 27^(1 - a;)V + 

7{ - 1 + 2a;(l ~ x){2 + x{l - x)} 

2(1 - 27 - 27^(1 - x)x + 27^(1 - x) V + 

7^(1 + 2a;(l - x){l + x{l - x)}]) 

—47a; 1 + X + 7^(3 — 2x)x — 

27^(1 - x)x^ - 7(1 + 4x2 - 2x^)" 



in the case of gluons. The x integral can easily be carried out but the result is simple only 
for the limiting values q = 0, /3 — > or /5 — >^ 1. The expressions for these cases are given in 
Sec. 4.2. 
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